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Abstract 

Since the experimental observation of quantum mechanical scattering phase 
shift in mesoscopic systems, several aspects of it has not yet been understood. 
The experimental observations has also accentuated many theoretical problems 
related to Friedel sum rule and negativity of partial density of states. We address 
these problems using the concepts of Argand diagram and Burgers circuit. We 
can prove the possibility of negative partial density of states in mesoscopic 
systems. Such a conclusive and general evidence cannot be given in one, two or 
three dimensions. We can show a general connection between phase drops and 
exactness of semi classical Friedel sum rule. We also show Argand diagram for 
a scattering matrix element can be of few classes based on their topology and 
all observations can be classified accordingly. 


1. Introduction 

A number of experimental m 1 n 10 n m and theoretical works [iini 
[TUmUinilTllIIllITlIITlITH] have studied the scattering phase shift in low 
dimensional mesoscopic systems where the electron dynamics is determined by 
quantum mechanics. Essentially, the asymptotic states are that of free particles 
in a quantum wire and the scatterer is embedded in the path of this quantum 
wire [Hiin]. Model specific approach so far has led to contradictory and con¬ 
fusing results [71 Uni EH EH EH HO] • Heiblum et al [7] quote that ‘There is by 
now vast theoretical evidence that the transmission phase depends on the spe¬ 
cific properties of the QD’s levels that participate in the transport’. However, 
in all these studies there is no reference to well established theorems and results 
on discontinuous phase changes and lapses known for a long time in the com¬ 
munity that studies classical wave trains miiiH. Such phase changes are still 
very important in mesoscopic physics as they are related to breakdown of parity 
effect 0 , interpretation of Friedel sum rule (FSR) [TOl EH EH ES] , relation to 
partial density of states (PDOS) |20l |2ll EH EH EH ES]; ^tc. that determine 
the thermodynamic properties of a mesoscopic system which is absent in case of 
classical waves. So in this work we first recapitulate why the scattering phase 
shift can be discontinuous as a result of Burgers circuit and then further use 
Burgers circuit to show its exact relation to density of states and partial density 
of states. Namely we can prove the reality of negative PDOS and prove the gen¬ 
eral regimes when semi classical FSR can be exact in a purely quantum regime. 
The advantage of using Burgers circuit is that these results can be shown to be 
very general and thus goes beyond all earlier works. 


Preprint submitted to Elsevier 


May 18, 2016 




Besides we will show that the scattering phase shift behaviour of different 
experimentally [TJ [H IH HI 0 HI [7] and theoretically [U HI [TOl HU [TTl [T31 EH 
EH EH EH EH studied mesoscopic systems can be understood from the Argand 
diagrams and analyticity. Each such system has its own peculiarities and so 
it is important to understand them with respect to a mathematical principle 
like Burgers circuit. Argand diagram for the scattering matrix element of these 
systems can be classified as, a) Argand diagram is closed, b) Argand diagram 
is open, c) Argand diagram encloses the phase singularity, d) Argand diagram 
does not enclose the phase singularity, e) Argand diagram is simply connected 
and f) Argand diagram is multiply connected due to the presence of sub-loops. 
Specific properties of the scatterer only matter to the extent that the Argand 
diagram changes from one of these to another. Among them some changes 
are topologically possible and others are not. Understanding how an Argand 
diagram changes from one to the other explains most of the puzzles. 

In section]^ we will analyse slips in the scattering phase shift [T1 ISl ITOl fT^ ■ 
In section we will show that FSR too can be understood from Argand dia¬ 
grams although its manifestation seems to be completely different for different 
potentials EH E3 El ESI EH- We will also show why FSR can become exact 
whenever there is a phase lapse. This is a physically counter-intuitive result 
that has been proven for particular potentials so far [H El El El Eg. We will 
show that this result depends on the properties of Argand diagram and hence 
very general and independent of the scattering potential. In section we will 
use Burgers circuit to prove the possibility of negative partial density of states 
in real mesoscopic systems. To show negative PDOS by explicit calculation 
for any particular realistic potential is virtually impossible and so never shown 
before. 

2. Burgers Circuit: an introduction 

Argand diagram is a plot of real versus imaginary parts of an analytic com¬ 
plex function and Burgers circuit is about phase changes and lapses being de¬ 
termined by phase singularities. Given a complex function t, if there is a phase 
singularity |22j in its complex plane, then one can specify the strength of the 
singularity as follows [77] . 


I = sgnim (Vt* x Vt) .n 


( 1 ) 


/ is a topological quantum number which is always conserved in an interaction. 
/ is a sign and so can be -1-1, —1 and 0. For a generalized ‘Burgers circuit’ (BC) 


dcj) = 27r/ 


( 2 ) 


where, (p = Arctan j ■ If the contour C does not enclose the phase singularity 
then / is 0. When the contour C enclosing a phase singularity is clockwise 
then / is -1 and when the contour C enclosing a phase singularity is counter¬ 
clockwise then / is -fl. A scattering matrix element is a complex function 
for which Argand diagram can be drawn and concept of BC can be applied. 
For the rest of the paper where we refer to scattering matrix element, we will 
mean the scattering matrix element of a quantum mechanical particle say an 


2 




Figure 1: Schematic Argand diagrams that exemplify Eq. §. (a) When the Argand diagram 
contour encloses the singular point at (0, 0) and (b) when the contour does not enclose the 
singular point. 


electron. Argand diagram for a scattering matrix element is generally counter¬ 
clockwise, as incident energy of the scattering particle increases. At the phase 
singularity, Re{t) = 0 and Im{t) = 0, implying the phase singularity is at the 
origin. Figs. dla) and Sb), shows schematic Argand diagrams for a complex 
function t. Fig. Sa), shows a typical counter-clockwise contour (ABDEA) of 
an Argand diagram enclosing the phase singularity at the origin. Hence for this 
case, I = -|-1. The contour trajectory (ABDEA) is concave throughout with 
respect to the singular point at the origin. Following Eq. ([^, the net change 
in phase in tracing ABDEA, in Fig. [^a), is 27r. Thus the phase monotonously 
increases for a concave trajectory. Fig. [^b) shows a typical contour (FGHJF) 
of an Argand diagram not enclosing the phase singularity at the origin. Hence 
for this case, 1 = 0. The contour has both concave (FGH) and convex (HJF) 
trajectories with respect to the singular point at the origin. Following Eq. ([^, 
the total change in phase in Fig. [^b) is zero. It is possible if the phase increases 
for the concave trajectory and decreases for the convex trajectory, the net phase 
change being zero. This can be easily verified by calculating (j) at any point 
{Re{t),Im{t)) of the trajectory and will be further demonstrated below. Thus 
for both Figs. Ba) and B^)) the phase change is determined by the singular 
point and the two follow the same principle. As a special case of the two, one 
can have a situation where the contour touches the singular point. In which 
case too the phase change can be understood from the same principle and will 
be explained later. Also if it is a closed contour as in Figs, [^a) and[2b), then 
Eq. ([^ is exact. In the real systems we will discuss in this work, the contours 
may not be closed. However, we will extend Eq. ([^ to understand such cases 
too. 

Whatever be the shape of a closed contour C, the phase change is given by 
Eq. As a consequence, the real and imaginary parts of t are not independent 
of each other but are related. Let us say the complex transmission amplitude 
be t{U) = where U can be any parameter like incident energy or 

gate voltage. Then 


Re {t{U)) 



U' -U 


dU' 


(3) 
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Im {t{U)) 



Re{t{U')) 
U' -U 


du' 


( 4 ) 


These are the well known Kramers-Kronig relations |30j . Another way in which 
the relation can be stated is in terms of Hilbert Transform E], 


ln\t{U)\ 



m') 

U' -u 


dU 


( 5 ) 


Ot{U) 



ln\t{U')\ 
U' -U 


dU 


( 6 ) 


Since / is a conserved quantity, by adding terms to a Hamiltonian (or details 
to the states in the scatterer) we cannot remove the phase singularities of the 
wave function in the complex plane. Phase changes are determined by the phase 
singularities. Depending on the interaction the Argand diagram can however 
change and so a theoretical understanding of the experimental data may not 
crucially depend on the details of the sample or model. Sample details can 
change the shape of the contour C, but as these theorems state, to understand 
the phase changes, we do not need all these details. Englman and Yahalom m 
had shown that the experimental data for scattering phase shift and scattering 
cross section of a quantum dot, are consistent with Hilbert transforms. We will 
show that the principles of analyticity and Eq. (H can be used to arrive at our 
results. 


3. Model Potentials 

In this section we analyse scattering phase shifts for different potentials that 
has been theoretically studied m so far, for analysing the experimental obser¬ 
vations [iiiiia. We intend to analyse w.r.t Eq. ([^ which has not been done in 
earlier works. Again as explained with Eq. Q, we do not need very complicated 
realistic potentials to understand the phase shifts but we need representative 
potentials that can be exactly solved and help us understand different aspects 
ofEq. ([|). 

3.1. Double delta function potential in one dimension 

Let us first consider scattering by a double delta function potential in one 
dimension (ID) schematically shown in Fig. [^a), that was studied in ref. m- 
Although a simple potential, it exhibits pronounced Breit Wigner (BW) reso¬ 
nances. We will use Eq. § to understand the scattering phase shift for this 
system and hence for BW resonances. The scattering potential for this system 
can be written as, 

Vi{x) = jiS{x) 

V 2 {x) = 72(5(a; - a) 

where 71 and 72 are the strengths of potential Vi and V 2 , respectively. The 
wave function in different regions marked /,// and III are (see Fig. [^a)), 

( e^^^ + re-^^=^ Jor X < 0, 
il}[x) = \ Ae*^^ + Be-^’^^Jor 0 < x < a, 

X > a. 
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Figure 2: (a) Schematic representation for scattering of electrons by a double delta function 
potential in one dimension. The direction of incident and scattered electrons are represented 
by arrows. The solid line represents a quantum wire with double delta function potentials 
at positions x = 0 and x = a respectively shown by cross (X) marks. 71 and 72 are the 
strengths of the potentials. The dashed lines represent the fact that the quantum wire is 
connected to electron reservoirs via leads, (b) Argand diagram for transmission amplitude for 
the double delta function potential, (c) Plot of transmission phase shift 9t versus ka and (d) 
plot of transmission coefficient |fp versus ka, for the double delta function potential using 
parameters e 7 ia = e 72 a = 40, a = 1, e = 1, 2me = 1 and h = 1. 
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Here r and t are the reflection and transmission amplitudes, k = 
is the wave vector and E is incident Fermi energy, t =\ t \ where, 

Of = Arctan is transmission phase shift and |f| = + Re{tY is 

transmission modulus. The Argand diagram for t is shown in Fig. [^b), where 
energy is varied to remain within the first Riemann surface. There is a phase 
singularity at the origin where t = 0. The Argand diagram encloses the singu¬ 
larity but is not closed in the first Riemann surface. Figs. [^c) and Hd) shows 
the transmission phase shift 9t and the transmission coefficient |tp, respectively, 
as a function of ka, using the same parameters as in Fig. [^b). 

In Fig. [^b), the contour starts from the origin where E = 0, goes first 
through point P and then through Q,R and S. The trajectory facing the 
singular point at the origin is concave throughout, and thus as discussed with 
Fig. 0 the phase increases continuously. This is evident in Fig. [^c) where the 
points P,Q,R and S are also shown at their respective values of ka. As the 
trajectory comes closer to the point of phase singularity, phase changes are very 
small and energy cost is very high. In Fig. [^b), the energy at the points marked 
P, Q, R and S are 8.614a^, 9.12a^, 23.6Ia^ and 37.58o^. Thus the energy change 
in going from P to Q (a large arc in the trajectory in Fig. [^b)) is very small, 
whereas the energy change in going from Q to P (a small arc in the trajectory 
in Fig. [^b)) is very high. The point R is very close to the singular point. Thus 
it costs a lot of energy as the Argand diagram trajectory tries to approach the 
point of phase singularity. 

3.2. Stub potential 

Another scattering potential often studied [13131] to understand the ex¬ 
periments [mai is known as the stub, which is an infinite one-dimensional 
quantum wire with a finite side branch. A schematic representation of this sys¬ 
tem is shown in Fig. [^a). This system is topologically not the same as a one 
dimensional system, as in this case the origin (0,0) (see Fig. [^a)) is connected 
to three other directions. This potential shows discontinuous phase drops by 
TT nni 131] similar to that observed in experiments BE]. Prior to that this 
discontinuous scattering phase shift was shown to cause breakdown of parity 
effect of single particle states [8]. However, if there is an energy scale associated 
with such a sharp phase change has remained a puzzle B- Heiblum et al B 
quote that ‘The appearance of a second energy scale in the phase jump between 
resonances also cannot be understood...’. They also mention that it suggests 
existence of an unusually large energy scale. We will use Eq. ([^ to understand 
why scattering phase shift can change discontinuously as some parameter is var¬ 
ied. Electrons are incident from left (see Fig. [^a)) with energy E. The thin 
lines represent one dimensional quantum wires with zero potential, while the 
bold line represents quantum wire with a finite potential V(y) given by, 

( 0,/or 0<y <li, 

V{y)= < iVoJor k <y<l, 

[oojor y>l. 

The potential V(y) is taken to be imaginary, as it allows us to make the Argand 
diagram trajectory approach and cross the point of phase singularity. This 
cannot be done with real potentials as / in Eq. ([^ is a conserved quantity. 


6 





Imaginary potentials are known as optical potentials [311 ESI ESI Ei] • They are 
often used to simulate the effect of decoherence and open systems. If the system 
changes from a closed one to an open one then I may be different for the two 
cases. When Vb is 0 then V{y) is 0 and the system is a closed conserved system. 
And when Vq is non-zero then V(jj) = iVo and it is an open system. The wave 
function in the different regions are given by, 


+re-^'^^,forx< 0 , 

te^^^,for X > 0, 


ip{x,y) = < 


+ Be-'^^yjor 0<y<li, 
Ce^Ay-h) + li<y <1, 


0,/or y = l. 


Here r and t are the reflection and transmission amplitudes, k = is the 

wave vector along thin lines, q = — iVo) is the wave vector along the 

bold line, and E is the Fermi energy. Solving the scattering problem using Grif¬ 
fiths boundary conditions [3111331 |5H |3S] , that the wave function is continuous 
and the currents are conserved at the junction at (0,0) in Fig. |^a), we get r 
and t as a function of energy, E. 


1 ik-qi -2ikl-t 

ik+qi ^ 

1 1 ni^—qi ^ — 2ikh 

^ "I" •^ik+qi ^ 

( 7 ) 

I-f 

( 8 ) 


where, qi = q [cot {q{l — Zi))]. Transmission phase shift is given by 6t = Arctan^^^. 
Fig. [3[b) show the Argand diagrams for the transmission amplitudes, given by 
Eq. for different values of eVoZ, where energy is varied to remain within 
the first Riemann surface. The thick solid line is for eV^l = 0, the dot-dashed 
line is for cVqI = 0.5 and the dashed line is for eV^l = —0.5. There is a phase 
singularity at t = 0 marked as S in Fig. [^b). Figs. |^c) and [^d) show the 
transmission phase shift 6t and transmission coefficient |fp, respectively, as a 
function of kl for different values of eVoZ, using the same parameters and same 
notations as in Fig. ib). 

In Fig. I^b), if we draw a circle (thin solid line through PABQP) around the 
phase singularity, then Eq. ([^ implies the total phase change along the contour 
is 27r. If the shape and size of the contour is altered, the net phase change 
along the contour remains the same provided the contour encircles the phase 
singularity. The phase change while going from Q to P along the arrow in the 
contour is less than tt as the phase change in going from R to A is tt. In the limit 
when the radius of the circle P, A, B, Q is tending to zero, P approaches A, Q 
approaches B and P, A, B, Q all coincide with the singular point S. In this limit 
going from Q to P would imply a discontinuous phase change of tt. Therefore 
a trajectory that tangentially touches the singular point S like the thick solid 
line for cVqI = 0 in Fig. ib), will exhibit a discontinuous phase change of tt 
at S. This phase change can be seen in the solid line for eVol = 0 in Fig. ic) 

(at point S marked at the same value of kl as in Fig. ib)), where we have 
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Figure 3: (a) Schematic representation for scattering of electrons by a stub potential. The 
direction of incident and scattered electrons are shown by arrows. The potential, represented 
by the bold line along y-axis, is V{y) = iVo- (b) Argand diagram for transmission amplitude 
for different values of eVb^ where, e is electronic charge and I is the length shown in (a). The 
thick solid line is Argand diagram for the case when eVol = 0, the dot-dashed line is that for 
eVol = 0.5 and the dashed line is for eVol = —0.5. Here li = .5/, / = 1, e = l,2me = 1 and 
h = 1. (c) Plot of transmission phase shift 9t and (d) plot of transmission coefficient as 
a function of dimensionless wave vector kl, taking the same notations and parameters as in 
Fig. I^b). 









































plotted the transmission phase shift 6t versus kl. We have already discussed 
in the previous section that as the Argand diagram trajectory approaches the 
point of phase singularity, the energy cost is very high. So, it is surprising 
that in a mesoscopic system, one can have such a high energy scale and one 
can realize discontinuous phase drops that are generally not seen for classical 
wave transport. In fact, it is not a real energy scale but an effective energy scale. 
Scattering by a stub of length I can be mapped into a problem of a delta function 
potential in one dimension where the strength of the delta function potential is 
kcot{kl). That is, {x) = kcot{kl)5{x) has same reflection amplitude r and 
transmission amplitude t as a stub of length 1. This is an effective potential that 
the electrons encounter while the real potential is Vq = 0. At kl = tt, a. small 
change in fc or Z means a very large change in the effective potential V'^^^{x). 
Such large effective energy scales are known in other areas of condensed matter 
physics, for example, the effective electron mass becomes oo at the band edge. 

For eVol = 0.5, the Argand diagram trajectory shown by the dot-dashed 
line in Fig. ib), intersects the Re{t) axis at point N at kl = n. The Argand 
diagram trajectory facing the phase singularity at S is concave throughout. So, 
Eq. ([^ implies a monotonously increasing phase. This phase behaviour can be 
seen in the dot-dashed line in Fig. [^c) where we have plotted the transmission 
phase shift 0* versus kl and marked N at kl = tt. For eVol = —0.5, the Argand 
diagram trajectory shown by the dashed line in Fig. ib), intersects the Re{t) 
axis at point M at kl = tt. The Argand diagram trajectory facing the phase 
singularity is partially convex and partially concave. As discussed with respect 
to Eq. Q , the phase increases for the concave part and decreases for the convex 
part of the trajectory. This phase behaviour can be seen in the dashed curve in 
Fig. i c) where we have plotted the transmission phase shift 9t versus kl and 
marked the point M similarly. Therefore, if eVgZ is continuously changed from 
0.5 to —0.5, then the point N moves to M crossing the point of phase singularity. 
When eVo^ becomes 0 we get the solid curves in Figs. ib),ic) and[^d) This 
implies that with an imaginary potential, just by changing sign of the potential 
it is possible to cross the singular point. / for dot-dashed line is 1 and that for 
the dashed line is 0. In other words by changing an imaginary potential we can 
make I change from 1 to 0. This is difficult with a real potential. If we tried 
to cross the singular point with real potentials it would cost an infinite amount 
of energy. We have thus discussed an effective real potential that can make 
an Argand diagram trajectory tangentially touch the point of phase singularity. 
Any effective potential that can make Argand diagram trajectory cross the point 
of phase singularity is not known. 

In Fig. I^c), we observe two different types of phase drops. For eVol = 0, 
we get a discontinuous phase drop (shown by the solid curve in Fig. [^c)) and 
for eVol = —0.5, we get a gradual phase drop (shown by the dashed curve in 
Fig. i c)). If we make eVol < —0.5, the point M in Fig. I^b) will shift more 
to the right and the phase drop will be less in magnitude and also less sharp. 
In Fig. [^d), using the same parameters and same notations as in Figs, [^b) 
and[^c), transmission coefficient is plotted as a function of kl for different 
values of eVol. The inset in Fig. ^d) shows the transmission coefficient 
in exponential scale in the region around kl = tt (shown by dotted circle). At 
kl = TT, the thick solid curve for eVol = 0 goes to zero, while the dot-dashed 
and dashed curves for eVol = 0.5 and eVol = —0.5 respectively, go through a 
non-zero minima and are very close. The phase behaviours for dot-dashed and 
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Figure 4: Schematic representation for scattering of electrons by a one dimensional delta 
function potential. The position of a delta function potential is shown by cross (X) mark. 


dashed curves are completely different as can be seen from Fig. [^c). 

Thus in this section, we have explained using Argand diagram and Burgers 
circuit, why phase (scattering phase shift) drops can occur discontinuously? 
Why such a phase drop can disappear or change from discontinuous to gradual? 
The gradual phase drop in the dashed curve of Fig. [^c) at M is related to the 
analytic property of complex transmission amplitude and how the trajectory 
encloses the singularity. Whether the gradual drop is sharp or slow, depends 
on the distance SM (Fig. |^b)) at which the convex trajectory intercepts the 
Re{t) axis. 


3.3. Single channel quantum wire 

The next scattering potential we consider, is a delta function potential in a 
single channel quantum wire [miss]. This essentially means, there is a single 
propagating channel while all other channels are evanescent. These evanescent 
channels are characteristic of quasi one dimension (QID), and make this scat¬ 
tering potential completely different from that of a delta function potential in 
one dimension [311 [Ml [3S]- This system gained relevance with respect to the 
experiments of Schuster et al. [T], Yang Ji et al. mm. etc., because it too 
shows phase drops [S] like those observed in case of the stub. The delta func¬ 
tion potential can create a bound state in the continuum and it is this bound 
state that non-trivially affects the scattering. Any other potential in quasi one 
dimension that can sustain a bound state will produce similar effects. Essen¬ 
tially, such a bound state cause a Fano resonance which is at the heart of the 
features observed for this potential. In this section we will explain the scattering 
phase shift for this potential from Argand diagram. The system is shown in Fig. 
[^ The quantum wire is shown by solid line with a delta function potential at 
position (0, yi) shown by a cross mark. W is width of the quantum wire. The 
dashed lines represent the fact that the quantum wire is connected to electron 
reservoirs via leads. Electrons are injected from the left reservoir into the left 
lead. The electrons are allowed to propagate along x direction, but confined 
along y direction. The confinement potential in the leads is taken to be hard 
wall potential and is given by. 


Vciy) 


oojor \y\ > 
0,/or |y| < 


The direction of propagation is shown by arrows. The scattering potential shown 
by cross mark is given by, 


V{x,y) = 'y6{x)d{y - yi) 
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Figure 5: (a) Argand diagram for transmission amplitude and (b) plot of transmission phase 
shift Otii versus EW^ for scattering by a delta function potential in a single channel quantum 
wire. Here e')W = —1.5, yi = 0.21VF, e = 1,H/ = 1 and we have considered 500 evanescent 
modes. 


Here 7 is the strength of the delta function potential. The asymptotic wave 
function in different regions are shown in Fig. One can solve the scattering 
problem [SS] to find, 


I'll = - 


'2ki 


1 + En>2 ^ +* 2 fc 7 


(9) 


tl 1 — 


1 + I2n>2 ^ 

1 + En>2 isr + 


( 10 ) 


Here rn and tn are the reflection and transmission amplitudes and r„m is given 

by, 

/ w\ 1 / w\ 1 

( 11 ) 



rZTT 

/ 1 F \1 


miT 

/ W \1 

7 sm 



sin 




where, m and n are integers. ki = ^~ is the wave vector for the prop¬ 
agating channel, denotes sum over evanescent modes, k„ = 
where n takes values 2, 3, ...00 and E is the incident Fermi energy. Transmission 
phase shift is given by, Otn = Arctan j ■ 

The Argand diagram for transmission amplitude tn is shown in Fig. (Ha). 
There is a phase singularity at the origin where, tn = 0 (shown by the point 
S). Fig. I^b) shows transmission phase shift as a function of energy, for 
the same parameters as in Fig. [^a). In Fig. [^a), the trajectory (shown by 
thick solid line) starts from the origin (point marked S) goes upto point Q, and 
traces back the same path to pass the origin making SQS a closed contour. The 
direction of the trajectory is therefore shown by a double headed arrow. At the 
energy, where the trajectory goes from Q to S' and touches the point of phase 
singularity at origin, i.e. point S, we expect a discontinuous phase drop of tt, 
following the same argument as in the case of the stub. This phase drop can be 
seen in Fig. gb), where the points Q and S are also marked at their respective 
energies. Thus the discontinuous phase drop is a natural consequence of Eq. 

mi¬ 


ll 

















Figure 6: Schematic representation of scattering of electrons by a three prong potential. The 
direction of incident and scattered electrons are shown by arrow heads. The potential is 
non-zero along the bold lines of lengths /, I 2 and I along —x^-\-y and +x axes, respectively. 


3.4- Three prong potential 

We now consider another potential called the three prong potential [TB] . This 
potential will help us to demonstrate other non-trivial aspects that follow from 
Eq. A schematic representation of the three prong potential is shown in 

Fig. 1^ The thin lines represent one dimensional quantum wires with potential 

V = 0, and the bold lines represent quantum wires with non zero potential, i.e., 

V ^ 0. The arms of the prong are labelled as 1, 2 and 3 as shown in Fig. 
The electrons are considered to be incident from left, the direction of incidence 
being shown by arrows. The wave function in the different regions are given by. 


ip{x,y,z) 


'e'^k{x+l) ^ ^ _i^ 

-I <x<0, 
Ce^’iy + De-^yyjor 0<y<l2, 
Jor 0 < x<l, 
t 2 iFk(y-i,) Jor y > l 2 , 
t^^e^Hx-i) jor X > 1. 


where fc = y is the wave vector along the thin lines, <7 = y {E — V) is 

the wave vector along the bold lines and E is the Fermi energy. Here rn is the re¬ 
flection amplitude for electrons incident from 1 and reflected back to 1 , ^21 is the 
transmission amplitude for electrons incident from 1 and transmitted to 2 and 
^31 is the transmission amplitude for electrons incident from 1 and transmitted 
to 3. These scattering matrix elements can be solved using Griffiths boundary 
conditions [SH [3311331 [33] . The respective transmission phase shifts are given 
by, = Arctan^^^, = Arctan-E^ and 0 * 3 , = Arctan-E^. 

Fig. I^a) shows the Argand diagram for transmission amplitude t 2 i- There 
is a phase singularity at the origin, where t 2 i = 0. In Fig. [^a) the trajectory 
of Argand diagram for t 2 i starts from the origin, goes through P and then 
through Q, R and S following counter-clockwise direction shown by arrows. The 
trajectory is concave throughout and Eq. implies monotonously increasing 
phase. This monotonously increasing phase can be seen in Fig. , where 
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Figure 7: (a) Argand diagram for transmission amplitude i 2 i and (b) plot of scattering phase 
shift 6 t 2 i S'S a function of kl varying the wave vector from 0 to 12. (c) Argand diagram for 
transmission amplitude fsi and (d) plot of scattering phase shift as a function of kl varying 
the wave vector from 0 to 12. For all the figures, I = 1, I 2 = 51, e = 1 and eVl = —1000. 


scattering phase shift 0 t 2 i is plotted as a function of kl and here also the points 
P, Q, R and S are marked at the corresponding values of kl. 

The Argand diagram for t^i shows something interesting. This Argand di¬ 
agram is shown in Fig. [^c). There is a phase singularity at the origin, where 
tsi = 0. In Fig. [^c), the trajectory of Argand diagram for ^31 starts from ori¬ 
gin, goes through P and then through D, A, B, Q, F, R and S, following counter¬ 
clockwise direction shown by arrows. Here interestingly, the trajectory develops 
a sub-loop ABQFA. This sub-loop results in a convex arc BQF in the trajec¬ 
tory that does not go through the origin. As explained earlier, there will be a 
gradual phase drop whenever such a convex arc is observed, following Eq. ([^. 
This can be seen in Fig. Bd), where scattering phase shift is plotted as a 
function of kl and here also the points P,A,Q,R and S are marked. Presence 
or absence of such a sub-loop has no consequence on the line integral of phase 
along PDABQFRS. This is because the contribution to the line integral com¬ 
ing from the sub-loop ABQFA is 0 and its presence or absence has no bearing 
on the value of I. So, such a sub-loop as ABQFA in Fig. [^c) can appear 
or disappear as some parameter is varied as will be demonstrated in the next 
section. 
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4. Injectance and Priedel sum rule 


Local partial density of states (LPDOS) is defined as [101 HSl 123 HZ] 


p'(a,r,/3) 


1 

47rz 


f -j- ^^aj3 \ 


(12) 


Here, =| Sap \ is the scattering matrix element for electrons incident 

from channel j3 and transmitted to channel a and stands for a functional 

derivative with respect to the local potential V{r). Time spent at r by an 
electron going from channel /3 to a is given by m 

T'(a,r,/3) = —^p'(a,r,/3) 

|Sa/3r 


Therefore, time spent by an electron going from channel /3 to a, within the 
scattering region is given by m 

T(a,/3) = — dr^p'{a,r,l3) (13) 

|^a/3| J screg 

where, ‘screg’ stands for scattering region. For a mesoscopic sample coupled to 
leads, this scattering region is by definition the sample m- So, partial density 
of states (PDOS) of a mesoscopic sample is defined as [27] 


p{a,/3) 
or, p(a,/3) 



(14) 

(15) 


PDOS are quite physical and manifests in a variety of experimental situations 
in mesoscopic systems [1311111Z]. For example, decoherence in the scattering 
region is proportional to the time electrons spend in the scattering region. As 
another example, consider a sinusoidal voltage of frequency ui, Vpito) applied at 
incident lead (3. The current measured at lead a will be m, 


la(uj) = Gap(uj)Vp(uj) 


(16) 


where, Gq,^(w) is the dynamical conductance matrix and is given by PZ] , 

Gap{<jj) = G^p — iioEap + + 0(a;^) (17) 

G", is the dc-conductance matrix. is proportional to io and governs the 

displacement currents and is given by [27] . 

Eap = e^p{a, l3)-e^ J dr'p{a, r') J drg{r, r')/o(r, /3) (18) 


where, g{r,r') is the effective interaction potential. All these experimental sit¬ 
uations explicitly involve /3 as the incoming channel and a as the outgoing 
channel. All analysis of p{a, P) will be made w.r.t the R.H.S of Eqs. (14) and 
(15). p{a,P) is theoretically undefined when the experiment does not clearly 
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involve an incoming channel and an outgoing channel. This has to do with the 
fact that quantum mechanics is necessary (to explain experimental observations) 
but not sufficient. There are many alternate approaches that give similar results 
as quantum mechanics and merits and demerits of such alternate approaches is 
a never ending discussion that we will avoid in this work. The fact remains that 
results deduced from quantum mechanics has never been violated in an experi¬ 
ment. So, whatever we prove for the mathematical expression on the R.H.S of 


Eqs. (14) and (15) will have consequences on experimental observations. 

The integration over r in Eq. ( [l4| can easily be done for a global change (for 
all r in the sample as well as in the leads) in V (r) by a constant amount e, i.e. 
SV (r) = e for all r. Such a constant global increase in potential is equivalent to 
decrease in incident energy E, i.e., 


and, therefore. 


5V{V) 


/ dr^-—— 

'sample 0 \/(r) 


d 

dE 


d 

dE 


(19) 


( 20 ) 


is expected to work in the semi-classical limit [5^. So, from Eq. (14) 


p{a,P) = 


_L f 


dsl 


Airi 


sample 


^al3 


,5E(r) 


Sa/3 


01/3 


SV{r) 


drHs^ 


Ss 


Ot^ 


5s\ 


4,711 


global 


(Using Eq. @) 


Q:/3 


or, p(a,/3) « — 




dSn 


Airi I dE 


dE 


On simplifying we get. 


p(«,/3) « ^ 




dE 


( 21 ) 


It is known that 


dOs 


dE 


^ can be negative. Concluding p(a, j3) to be negative 
when R.H.S of Eq. (21) is negative is completely wrong. Eq. is an approx¬ 
imate equality, which implies if R.H.S is negative, the L.H.S is not necessarily 
negative. Although PDOS p{a, jd) can be in principle negative as can be demon¬ 
strated from R.H.S of Eq. (14) for some strictly ID simple potentials (like a 
square well and square barrier) in very low energy regime. ID potentials are 
an idealization and not physical and neither one can go to the necessary low 
energy regime in an experiment. This work is based on the realization that 


some general conclusions can be drawn about the R.H.S of Eqs. (14) and (15) 


using the properties of Argand diagram and Burgers circuit which is true for 
any potential in ID or in QID. 

The problem of negative PDOS was theoretically studied in two different 
ways strictly in ID. One is negative PDOS and the other is negative time scales. 
Measured time scales in real systems did not lead to any unique physical un¬ 
derstanding and a review on the topic can be seen in ref. |30]- Ref. [JD] did 
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conclude that p{a,P) or T{a,l3) are physical and later on found to be so P7] . 
But regimes where they become negative has not received any theoretical at¬ 
tention beyond ID. In other words negative p{a,j3) has not yet been shown for 
any physical system. 


Although PDOS as defined in Eq. (14) can be negative, they add up to give 


the correct DOS m which is positive. One can sum the PDOS (given by Eq. 
14) over a to get injectance PHI HZl US] of lead (3, 




dr'^ 


6s 


' sample 


^oi(3 


Oi^ 


5sl 


5V{v) 


Sal3 


OL^ 


,5E(r) 


( 22 ) 


This is a measure of current delivered by the electrons incident along lead /3 
and outgoing through all the leads. Similarly injectance can be defined for all 
possible leads and they are completely independent of each other. Eor example, 
the scattering problem depicted in Fig. |^can define p(l). To get p{2) one has 
to solve a completely different scattering problem where the incident particle is 


from lead 2. Using Eq. (20) we get the semi-classical limit of Eq. 




ds. 


OL^ 


ds 


,t ' 

a/3 


Sa/3 


(23) 


Summing p(/3) over the /3 independent channels, we can obtain density of states 
(DOS) p{E), i.e., from Eq. (@, 






OL 0 


> sample 


^a/3 


(5U(r) 


Sap 


and in the semi-classical limit given by Eq. (20), we get 

dsap 


1 




ot 


Further simplification of R.H.S in Eq. (25) gives 


ds 


dE dE 


dU(r) 


(24) 


a/3 


(25) 




(26) 


This is Friedel sum rule (FSR), where 0f{E) = :^-J^log{det[S]) is the Friedel 
phase, S is the scattering matrix and p{E) = is density of states. Since 

injectance of all leads are independent while they add up to give DOS, it is 
important to understand injectance in order to understand FSR. So we will 
restrict our study to injectance. 

The potentials in sections [3.2[|3.3| and |3.4| are typical examples of mesoscopic 
systems, and as reported in ref. [iniiiiiiisi, FSR (or injectance) manifests in 
these systems, in different ways. Semi-classical regime being expressed by Eq. 


(20) does not seem to be sufficient. Sometimes, FSR (or injectance) is exact at 


all energies (for example the stub |10j l and sometimes it is exact at an energy 
where quantum fluctuations dominate m- There is a huge amount of system to 
system variation. It has been proved very generally that when the phase drops 
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of TT are discontinuous like that in the solid curve in Fig. [^c), then at the energy 
corresponding to this drop, semi-classical injectance will become exact [15) . But 
when the phase drops are gradual like the dashed curve in FigJ^c), this has not 
been proved in general but shown for particular cases [El 129)7^ Burgers circuit 
will help us to derive in this work the general connection between gradual phase 
drops and exactness of semi-classical injectance for any arbitrary potential that 
can exhibit such a gradual phase drop. 

We will first prove that in realistic mesoscopic systems one can definitely 


observe negative PDOS. We will show that the R.H.S of Eqs. (14) and (15) can 


become negative for a realistic mesoscopic system. We will also show that such 
a conclusion cannot be drawn in ID that has been extensively studied before 
m- Also we will show that when there are such negative slopes in scattering 
phase shift of mesoscopic systems then semi-classical FSR can become exact 
in a quantum regime. Once again our proofs will depend on Argand diagram 
and Eq. (|^ and so our proof will be general and not depend on the specific 
properties of the scattering potential. 

For this we first make the connection between Eq. ([^ and injectance. When 
S matrix is 2 x 2 (for example, the cases of double delta function potential 
(Fig. I^a)), stub (Fig. [^a)), single channel quantum wire (Fig. |^), Eq. ( |2^ 
simplifies to. 




Trp{E) 


(27) 


Here t is the transmission amplitude and 6t = Arctan ■ It has been proved 
that this equation is valid even if 9t is discontinuous with E because 9t is analytic 
[28] wherein the R.H.S derivative of 9t is the same as the L.H.S derivative at the 
discontinuity. Suppose when the energy is varied from 0 to i?i, then the Argand 
diagram for a typical scattering matrix element t traces a closed contour C. For 
such a case one can state that, Eq. ([^ takes the form 


d9t = 2 'kN{Ei) 


(28) 


where, N{Ei) is number of states (obtained by integrating DOS p{E) from 0 
to El) below energy Ei and is to be identified with the conserved quantity 
I. Comparing with Eq. ([^ we see (j) = 9t and I = N{Ei). If C happens 


to be a completely closed contour, Eq. (28) is exact as it is equivalent to Eq. 
(H). Any complex function or its phase has to satisfy Eq. <§ and a scattering 
matrix element is no exception provided it is analytic. This analyticity is the 
basic requirement for Eq. ([T^ (PDOS) or Eq. (22) (injectance). However, 


in most cases of scattering problems, C is not completely closed and when C 
is not completely closed one cannot expect any conserved quantity. In Fig. 
Sb), there is a phase singularity at the origin and the contour enclosing the 
singularity is not closed in the first Riemann surface. One can restrict the 
discussion to the first Riemann surface to understand the injectance. When the 
contour continue to the second Riemann surface, then the contour integral starts 
including contribution from the second phase singularity in the second Riemann 
surface. And then one has to extend the discussions here to include the effect of 
the second singularity. This does not change the arguments given here except 
that sometimes the error from the first Riemann surface can cancel the error 
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from the second Riemann surface which need not be a systematic behaviour. 
The contour C in Fig. [^b) starts from origin with zero energy and ends in the 
first Riemann surface at point marked as S, where the energy is Ei (say). It is 
now known that (see Eq. (6) in ref. |13)b 


II 

pEi - 

JC 

Jo 




dE 


dE 


dE 


(29) 


We will replace C by C', when the contour is not completely closed. G"^ is the 
advanced Greens function and is self energy due to coupling the system with 
the leads. One can then state that ImTr dE is the correction term 

for Eq. (28) when contour C' is not closed. This statement can be alternately 
justified as follows. When the self energy is independent of incident energy, then 
the contour C is closed as well as the correction term is zero implying Eq. (29) 
becomes Eq. (28). For a double delta function potential in one dimension the 
correction term is very important to consider. There, energy dependence of self 
energy can be seen very easily in the broadening of consecutive resonance peaks 
(shown in Fig. [^d)). Therefore, one can refine the statement in Eq. (20) to 
state, 


dE / 

C sample 




(30) 


where, we generate a closed contour C in the Argand diagram by varying energy 
E from 0 to El. Without this equality one cannot get Eq. (28) as can be easily 


verified. Now from Eq. (19) one can state. 


dE / dr^—— = f dE 

C J sample dV[T) Jq 


/ dr^- —— 
global dV (r) 


(31) 


In case of the solid line in Fig. ib), the potential everywhere is real and is 
0. It traces a closed contour in the first Riemann surface. Thus for this system 


Eq. (28) will become applicable. Explicit calculations of density of states m 
for the stub show this and so everything is consistent with Eq. ([^. However, 
for scattering by a delta function potential (Fig. in a single channel quantum 
wire we get a counter intuitive result. In this case, the contour of the Argand 
diagram (shown in Fig. [^a)) is closed in a special way and so we expect Eq. (28) 
to be applicable. But explicit calculations of density of states m show, that 
is not the case. It has been shown earlier that for this system the fundamental 
theorem of Biittiker-Thomas-Pretre (BTP) also breaks down [5T] due to the 
non-analyticity of scattering matrix elements. This is a consequence of the fact 
that delta function potential in QID incorporates a log divergence in scattering 
matrix elements. Hence this is a situation where Eq. (§ cannot be applied. 

For the three prong potential, shown in Fig. [6} t he scattering matrix is 3 x 3 
and the correct form of FSR is given by Eq. p6|. Whenever the scattering 


matrix has a rank greater than 2, the connection between FSR in Eq. (26) and 
Eq. ([^ is not straight forward. However, we can make this connection for each 
partial density of states (PDOS) and is shown below. As an example, let us 
consider the Argand diagrams for the three prong potential shown in Figs. Ha) 
and Be). None of the Argand diagrams (e.g. Figs. c)) are closed in the 

first Riemann surface. Let us, for example, consider the Argand diagram for 
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Figure 8: ^ ^ function of kl for the three prong potential. Here / = 1, ^2 = 5/, e = 1 

and eVl = —1000. 


tracing a contour C' {PDABQFARS shown in Fig. [^c)) as energy is varied 
from 0 to ifi- We can show that (see Appendix A), 


d9u 


2 . r 

Jo Fail 


(32) 


Now we can again state that for any closed curve like ABQFA in Fig. He), Eq. 
(32) is exact. That is, 


dOu 


ABQFA 


= 27r ® 

JABQFA 


d3,l)(A) 


dE 


(33) 


This statement can also be alternately justified as follows. For a closed contour, 
L.H.S of Eq. (pSl) is zero. Also for a closed contour R.H.S of Eq. (|33|) will be 


zero as shown in Appendix B for any p(a, fd) inside a closed contour integral. 
Note that in Appendix B the expression that we have used for p[a, fd) is given 
by R.H.S of Eq. (15). Hence it follows that R.H.S of Eq. (33) is also zero, where 


p(3,1) or any p{a,/d) is given by the R.H.S of Eq. (15). The arguments below 
although stated for p(3,1) is therefore true for any p{a,ld){E) inside a closed 
contour integral. Now, 


ABQFA 




de, 


^-JJ-dE = 0 


ABQFA 


dE 


(34) 


ddt. 


In Fig. H is shown in the energy range covering the sub-loop ABQEA 
of Fig. ph). As is implied by Eq. 


(34 


somewhere negative to ensure the area unc 


8 ) is zero. Similarly the R.H.S of Eq. (33) is zero implies that will 


d9t. 


Jff- is somewhere positive and 
er the curve (shaded region in Fig. 


so be positive as well as negative in certain energy values (or kl values). Thus 


p(3,1) as given by the R.H.S of Eq. (15) is conclusively negative in some energy 
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values. Although we have considered the there prong potential as an example 
the proof is valid for any potential whose Argand diagram shows a continuous 
phase drop due to a sub-loop as it is due to the fact that the sub-loop traces a 
closed contour. 

We will now show how negative slopes in scattering phase shift of mesoscopic 
systems are fundamentally different from that studied earlier [ID] in ID. Note 
the negative slope at point P in Fig. m- This kind of negative slope at very 
low energies can arise for scattering in ID and one can easily check this for a 
square well potential. In terms of our analysis we understand the negative slope 
at point P in Fig. [Tjjd) due to a convex trajectory at P in Fig. I^c) which 
is originating due to the fact that the Argand diagram starts from the origin 
and behaves anomalously as the trajectory starting from the origin is neither 
clockwise nor anti-clockwise with respect to the singular point (i.e. origin). 
See the expanded Argand diagram trajectory shown in the inset of Fig. [^c). 
The trajectory moves up, turns around and moves down to become convex in a 
small energy window and then winds around the origin anti-clockwise. Although 
is negative at P in Fig. [^d), there is no conclusive evidence that PDOS 
p{3, 1) is negative at energy corresponding to point P. Such a negative slope is 
fundamentally different from the negative slope at Q in Fig. [^d) that originate 
from a closed sub-loop ABQFA in Fig. [^c), that we encounter only in QID 
and mesoscopic scattering. We have shown that the sub-loop seen in Fig. [^c) 
implies the presence of this negative slope in scattering phase shift 9t^-^ and also 
implies PDOS p{3, 1) is negative. We have shown that an Argand diagram for 
such a scattering matrix element curls around and forms a sub-loop without 
violating the topological constraints of Eq. (H- The line integration along the 
sub-loop ABQFA in Fig. [^c) does not contribute to the line integration over 
the trajectory PDABQFARS of Fig. [^c) or I in Eq. ([^ is unaffected by the 
presence or absence of a sub-loop in the closed contour C. Negative slopes of 
this second type (i.e. observed at Q in Fig. m) that we have discussed here 
in-fact can appear or disappear very easily and can be found at much higher 
energies. In Fig. we have shown the Argand diagram for Qi upto very 
high value of energy (or, kl). We can see many sub-loops which again implies 
the presence of negative slopes in scattering phase shift and negative PDOS. 
Sometimes there is a cusp in the Argand diagram and such a cusp means a 
sub-loop has disappeared [21]. Thus the phase drop will also disappear as we go 
to such energies. Disappearance of a sub-loop can be demonstrated by varying 
any other parameter like V,l,l 2 of the three prong geometry. 

p(3,1) being negative is a counter intuitive feature of quantum mechanics 
and can have interesting physical significance. Obviously, ac-response of a meso¬ 
scopic system will change drastically if p{a, P) in Eq. ( [I^ changes sign. Also, it 
means an electron that is incident along lead 1 and transmitted to lead 3, dwells 
in some negative number of states (PDOS is negative) inside the scatterer. To¬ 
tal charge being electronic charge times number of states will be positive for 
these negatively charged electrons. So other electrons that are incident along 
lead 1 and transmitted to lead 2 or reflected back to lead 1 will be attracted by 
this positively behaving charge of electrons going from 1 to 3. This could be 
the explanation for the electron-electron attraction observed in the numerical 
experiment of ref. |42] , where no explanation could be given. In the next section 
we will argue that this could have been also observed in an experiment. 
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Figure 9: Argand diagram for transmission amplitude tai, varying kl from 0 to 50, for the 
three prong potential. Here / = 1, ^2 = 5/, e = 1 and eVl = —1000. 


Now let us try to understand if there is a general connection between negative 
slope and injectance becoming exact as observed in some earlier works [29j for 


specific potentials. Eq. (33) holds for the integrals and does not imply equality 


of the integrands. However, using Eq. (211 we can write 


p(3,l) s 

1 U |2'^^t3l 

27r^ dE 

(35) 

p(2,l) « 

27r^ dE 

(36) 

p(l,l) R 

1 

■ 2,1"“' dE 

(37) 


Finally using Eqs. (35), (36) and (37), we can write 

p(l,l) + p(2,l)+p(3,l) 


1 

27T 


rii 


2 dOrii 

dE 
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Ih 


>d9t. 


dE 


27r^^^^ dE 


(38) 


L.H.S of Eq. (38), i.e., X)„p(a,/3) is well known as injectance as defined by 


the R.H.S in Eq. (22). The R.H.S of Eq. (38) is the semi-classical limit of 


injectance. The topological interpretation of Eqs. (35), (36), (37) in terms of 


Argand diagrams leading to Eq. (38) is very useful. The correction term to Eq. 
((38) is known [511135] i.e., 


p(l,l)+p(2,l)+p(3,l) = 


1 

27r 


^ dE 


,d9, 


'-211 


*21 


dE 




de,., ^ 


dE 


hk^ 


(39) 


Let us say that, 
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Figure 10: (a) Argand diagram for reflection amplitude rn and (b) plot of reflection phase 
shift 9rii , versus kl for the three prong potential. Here / = 1, /2 = 5/, e = 1 and eVl = —10000. 



Figure 11: Plot of exact injectance p(l)^ (solid line) and semi-classical injectance p(l)'^ (dashed 
line) as a function of kl for the three prong potential. The peaks in the injectance are shown 
separately, (a) shows the first peak, (b) shows the second peak, for the same parameters as in 
Fig. |10[ The insets show the magnified curves at points M and N. 


p(l)^ = p(l,l) + p(2,l)+p(3,l) 


(40) 


is exact injectance, defined by the R.H.S in Eq. (22). And 
1 I ,cyd0r^^ I. 10^0^21 




rii 
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dE 


\t2 


dE 


Ihi 
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dE 


(41) 


is generally referred to as semi-classical injectance. Eq. (39) implies that 
p(l)® and p(l)® will be equal at energies, where the correction term is zero 
or \rii\sin{6rii) = 0. According to the arguments of Leavens and Aers [23], 
in the semi-classical limit |rii| —>■ 0 and p(l)® = /^(l)®. But in the cases of 
studies on injectance or FSR in mesoscopic systems [151129], I rill 7 ^ 0 at the 
energy where p(l)® = p(l)®, but not very consistently as already discussed in 
the paragraph after Eq. ([2^. We will show how the Argand diagram topology 


is responsible for this behaviour and therefore provides a general understand¬ 
ing. In Figs. [^a) and Mh), we have shown the Argand diagram and phase 
shift for scattering matrix element rn, respectively, kl is varied from 0 to 12 
in both the plots. The Argand diagram in Fig. [^a) is restricted to one side 
of the phase singularity (i.e. origin) and in the first Riemann surface resulting 
in sub-loops. This will naturally mean that, the contour has both concave and 
convex parts in the trajectory. Scattering phase shift in Fig. HKb) increases 
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Figure 12: Plot of reflection phase shift flrn versus kl for the three prong potential. Here 
Z = 1, Z 2 = 51, e = 1 and eVl = —1000. 



Figure 13: Plot of exact injectance p(l)'^ (solid line) and semi-classical injectance p(l)^ (dashed 
line) as a function of kl for the three prong potential. The peaks in the injectance are shown 
separately, (a) shows the first peak, (b) shows the second peak, for the same parameters as in 
Fig. |12| The insets show the magnified curves at points F and G. 


with kl, reaches a peak value and then drops to become tt at M. The pat¬ 
tern repeats as kl increases and the scattering phase shift becomes tt again at 
N. Therefore the correction term \rii\sin{0r^^) to semi-classical injectance (Eq. 
(39)) is zero at M and N. Thus at M and N, the semi classical injectance (Eq. 
(41)) is exact. The exactness of semi-classical injectance is shown in Fig. [TI](for 
clarity see the inset) at points M and N corresponding to same kl values as in 
Fig. |10[ For monotonously increasing phase, Argand diagram extends to higher 
Riemann surfaces and line integrals include the effect of other singularities in 
higher Riemann surface. Of course phase can be integral multiples of tt (i.e., 
27r, Stt, ..), but for open Argand diagram trajectories as argued before, this is 
not a consistent behaviour, except in some simple one dimensional scattering 
problems. This inconsistent behaviour can be for example, checked for an one 
dimensional Aharanov-Bohm ring with different arm lengths. Drops in 0^1 re¬ 
sulting in 0^11 being tt and sm(0m) = 0, leading to semi classical injectance 
being exact can be understood from the Argand diagram in a single Riemann 
surface, involving a single phase singularity in the line integral and hence is not 
an accident. Drops in 0^1 1 is a pure quantum mechanical behaviour and hence 
exactness of semi classical injectance at the energies corresponding to the phase 
drops, is counter-intuitive. As argued before, these drops coming from sub-loops 
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Figure 14: Transmission coefficient and scattering phase shift evolution in a QD as the coupling 
between the dot and leads is changed. This figure is taken with permission, from ref. [2]. 


in Argand diagram are tunable and can be removed by varying some parameter. 
In Fig. we plot the phase behaviour for the same system, after decreasing 
the potential V of the same system as in Figs. 10 and 11 The phase drop at 
point F as usual decreases to tt, and hence semi classical injectance is exact at 
this point which can be seen in Fig. (Ha). Point F is marked at the same value 
of kl as in Fig. But, in Fig. the phase drop at G is now not sharp enough 
to decrease to tt. It is due to tuning the potential, that the sub-loop has now 
reduced in area, and consequently the drop has also reduced and will eventually 
disappear as V is decreased further. At such point p(l)® is not equal to p(l)®, 
i.e. semi classical injectance is not exact in spite of a phase drop. This can be 
seen in Fig. [T^b)) where point G is marked at the same value of kl, as in Fig. 
|12| One would have expected that when we decrease the potential and make 
it weaker, semi-classical behaviour will be favoured. But on the contrary, for 
stronger potential the semi-classical injectance is exact at point N in Fig. [I^b), 
while for a weaker potential in the same system the semi-classical injectance is 
not exact at point G in Fig. [T^b). Therefore, the drops in scattering phase shift 
of mesoscopic system, originating from sub-loops in Argand diagram involves 
completely new physics. One has to discard the usual concept of semi-classical 
regimes wherein the de-Broglie wavelength of the electron is much smaller than 
the scale of the potential, mathematically expressed by Eq. (20) [20] . 


5. Comment on Experimental Observations 

We will try to construct the Argand diagram from the experimental data 
given in refs. mill, to see if something can be understood. Yang Ji et al. m 
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Figure 15: Argand diagram for the transmission amplitude, obtained from Fig. |14| In this 
figure the labelling A,B,C,D follows that in Fig. 


measured the transmission coefficient and transmission phase shift of a quantum 
dot (QD). The experimentally observed transmission coefficient and transmis¬ 
sion phase shift as a function of gate voltage, for different coupling strengths 
between QD and the leads, are shown in Fig. A phase drop is observed in 
Fig. [153 when the coupling strength is small. As coupling gradually increases, 
the dot enters the strongly interacting regime and the phase drop decreases and 
finally vanishes. This gradual phase disappearance follows the sequence D to A 
in Fig. We have seen in our theoretical calculations that such phase drops 
can disappear and there are two possible reasons. One possibility is that when 
the coupling to leads is changed then it introduces some de-coherence for which 
the Argand diagram trajectory crosses the singularity and / changes from 0 to 
1. The second possibility is that I remains conserved and the Argand diagram 
trajectory develops a sub-loop. Argand diagrams can be constructed from the 
experimental data given in Fig. and is shown in Fig. There is a lot of 
fluctuations in the data probably due to experimental error. However, one can 
roughly see that there is a sub-loop abca in Fig. that gradually decreases 
in area as we go from D to A and clearly disappears in A with the formation 
of a cusp. The experimentalists had thought that as coupling increases, Kondo 
effect sets in and the disappearance of phase drop could be due to that. How¬ 
ever, the phase behaviour seen in Fig. is not typical of Kondo resonance. 
Typical phase behaviour of Kondo resonance is that, it increases by ^ and forms 
a plateau. Such a behaviour was observed by the same experimental group [7] 
subsequently, in which they restricted the dot occupancy to one electron, and 
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Figure 16: (a) Typical transmission coefficient and scattering phase shift of a QD in the 

Kondo regime. This figure is taken with permission from ref. [?]• (b) Argand diagram 

for transmission amplitude constructed from the experimental data, (c) We hope future 
experiments will give such Argand diagrams in greater resolution and how they evolve. For 
example an Argand diagram that encircles the phase singularity (solid line) can develop lobes 
(dotted line) and finally reduce to the Argand diagram in (b) as Kondo effect sets in. 
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is shown in Fig. |16[ a). One may conclude that the disappearing phase be¬ 
haviour in Fig. |14| 4 has nothing to do with Kondo effect. We have already 
discussed in section that when Argand diagram of transmission amplitude 
shows a sub-loop, PDOS associated with the transmitting electrons is negative. 
The transmitted particles and reflected particles can attract each other and that 
can result in a bound state. This bound state is expected to be effective only 
in the strongly interacting regime, i.e. for the case of Fig. In the weakly 

interacting regime, we may not see any effect of this interaction induced bound 
state and the phase drops occur very generally in scattering by a mesoscopic 
system. We hope future experiments will focus on how Argand diagram evolves 
as one observes changes in scattering phase shift. 

6. Conclusions 

Argand diagram of scattering matrix elements are drawn for different model 
potentials and for a few experimental data. Several conclusions can be drawn 
from the topology of the Argand diagram without referring to the Hamiltonian 
or to the scattering potential. In ID, 2D and 3D, the Argand diagram trajec¬ 
tory encircles the phase singularity. But in mesoscopic systems we find that the 
Argand diagram develops sub-loops. The sub-loop does not enclose the phase 
singularity at the origin and hence topologically allowed. The sub-loop can 
therefore appear or disappear on varying some parameter. It does not matter 
what parameter (say E) is varied to obtain the Argand diagram or what pa¬ 
rameter (say V) is varied to make the sub-loop disappear. Many unexplained 
features so far can be explained by the appearance and disappearance of such 
a sub-loop. When the sub-loop appears there will be a gradual drop in the 
scattering phase shift and when the sub-loop disappear the drop will also dis¬ 
appear. Hence appearance and disappearance of phase drop is also very natural 
and poses no conceptual problem. Just as the sub-loop appears on varying some 
parameter, it can also grow in size as the parameter is varied. As the sub-loop 
becomes large and comes closer to the origin the phase drops also become large 
and sharp making the scattering phase shift decrease to tt and then injectance 
or Friedel sum rule becomes exact. This is very counter intuitive as the strong 
phase drop signifies onset of pure quantum mechanical behaviour while Friedel 
sum rule is expected to become exact in semi-classical regimes. Also we prove 
that whenever there is a sub-loop (big or small) there will be negative partial 
density of states. For example if there is a sub-loop in the range AA then there 
is also negative partial density of states in the range AA. Conclusive evidence of 
negative partial density of sates in a real system has never been reported before. 
Since all these results are drawn from the properties of the Argand diagram, 
the results are general and independent of the Hamiltonian or the scattering 
potential. The physics originating from sub-loops is completely new and upsets 
our way of understanding semi-classical behaviour. 
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Appendix A 


From Eq. (15), we can write 

MM) = 4 


Jsample Ok (rj 


In Eq. (A.l I using the well known semi classical approximation dv^ , 5 y (r) — 

1 12'^^*31 


[20], we get, (see Eq. (21)) 


M3.1) » 

Now, for a contour C" traced when energy is varied from 0 to Ei, 

^ dOt 


ic 


d-dt^i — 


7*31 

dE 


dE 


rEi 12 ^^*31 

27rr31l 


Using Eq. (A.2) 


de 


*31 


27r 


P(3,l) 

\t31? 


dE 


dE 


(A.2) 


(A.3) 


Appendix B 


We know that in a scattering problem increasing incident energy by dE is 
equivalent to decreasing the potential globally by a constant amount Ae, such 
that dE = —eAe, where e is particle charge that we will set to 1 to simplify our 
arguments. That is, the new potential is U'(r) = U(r) — Ae. Hence if we can 
generate a closed sub-loop in the Argand diagram by varying E, then we can 
also do so by globally changing the potential and for such a closed contour like 
ABQEA in Fig. [7]:c), 


i.e. 


ABQFA 

56. 


50s^, = 0 


ABQFA J global 


(5U(r) 


Aedr"^ = 0 


(B.l) 


Now we replace the global integration over r by an integration over the sample 
or the scattering region only, since we have seen that it can be done in case of 
closed contours or inside an integration of the type in Eq. (B.l). This has 
already been discussed in Eq. (31). 


Sa/3| 


6es 


' I I ® 

C dsample ol/(rj 


or, 
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rAe 


,2 


C J sample 

or, 27r 


(5U(r) 

p{a,l3) 

C l®Q;/3p 


= 0 
Ae = 0 


Using Eq. (15) 


Therefore, the R.H.S of Eq. (33) is justified if an electronic charge is multiplied 
to the numerator. 
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